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Exterior calculus is a powerful tool to search for solutions to the electromagnetic field equations, 
whose strength can be better appreciated when applied to work out non-trivial configurations. Here 
we show how to exploit this machinery to obtain the electromagnetic TM and TE modes in hollow 
cylindrical waveguides. The proper use of exterior calculus and Lorentz boosts will straightforwardly 

(~~) ' lead to such solutions and the respective power transmitted along the waveguide. 
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(N ; I. INTRODUCTION 

c : 

Waveguides are an excellent arena to practice relativity and calculus with differential forms, since these tools greatly 
help the understanding of electromagnetic problems and their solutions. Although most of textbooks teaching exterior 
calculus to physicists display some basic applications to electromagnetismji~— any of them stresses the power of this 
language by applying it to solve Maxwell equations in waveguides and cavities. This is rather disappointing, because 
the skill in applying both tools -relativity and exterior calculus- allows a clean presentation of the subject and an 

O y easy way for working out the field configurations, so highlighting the power of the language not only as a theoretical 
weapon but as a practical tool to solve problems of technical relevance. On the contrary, normally the vector language 
is used to display the field configurations in waveguides (however see Ref. [H), which is hardly a natural language for 
electromagnetism. Because of this reason, the vector approach tends to be rather tedious. Frequently, the lack of 
a proper geometric language reduces the explanation to a simple case like the waveguide of rectangular or circular 
■ section (however see Ref. ly). We will work out the problem of propagating waves in hollow cylindrical waveguides of 
arbitrary section by employing exterior calculus. We will start with TM (TE) (non-propagating) stationary modes, 
whose field structure is very simple: it is just an electric (magnetic) field along the waveguide and a magnetic (electric) 
field in the waveguide section. Then, we will turn these stationary solutions into propagating solutions by means of 
a Lorentz boost along the waveguide. Finally we will compute the transmitted power and emphasize its relativistic 
relation with the energy per unit of length in the waveguide. The calculus involving exterior derivatives, Hodge 
dualities and the generalized Stokes theorem will provide a straightforward way for building the solutions, since all 
£C) . the vector equations of the usual approach condense in just one equation in geometric language, so showing how 

^ ' powerful this machinery is. 
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Exterior calculus is the natural language for electrodynamics. 7,8 But not only electrodynamics greatly benefits 
from the compactness and simplicity of this geometric language. Also the developments in Hamiltonian mechanics*^ 
thermodynamics^ 1 '- Yang-Mills fields /^ 11 geometric (Berry) phases in quantum mechanics^ topological quantum 
fields as the Chern-Simons theory ; 1 ^ 11 gravity?^ symplectic geometry^ connections in fiber bundles^ etc. gain in 
clarity and depth when expressed through the tools of exterior calculus. There is a huge list of textbooks to learn 
exterior calculus. The reader is referred to Ref. li-4 for a first approach to the subject. Here we will briefly review 
the main operational features of the exterior derivative d and the wedge product A between differential forms. 

Any linear combination of coordinate differentials at each point of space is a 1-form field (whatever the coordinates 
are, Cartesian or not; even if the geometry is non-Euclidean). For instance, we may call r\ the 1-form 

rj = 3 x 2 y 7 dx + 5 y dz (1) 

If the space is 3-dimensional and (x, y, z) are the chosen coordinates, we say that {dx, dy, dz} is a coordinate basis for 
1-forms. The components of the above defined 1-form rj in this basis are r\ x = 3x 2 y 7 , r\ y — 0, r/ z — 5y. Generically, a 
1-form in a n-dimensional space is 

a = dx^ 1 (2) 

(Einstein convention is used). The superindex in dx^ 1 labels the n 1-forms of the coordinate basis, whereas the subindex 
in a.n labels the components of the 1-form a, which are functions of the coordinates. 
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II. A BRIEF REVIEW OF EXTERIOR CALCULUS 
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1-forms can be introduced in a geometric way as linear real valued functions on the tangent vector space to a 
manifold: they are covectors or covariant vectors*^ However, here we are not interested in the action of forms on 
vectors. Instead we will operate within the set of p-forms, which are defined as totally antisymmetric covariant tensors 
of p indexes (so p ^ n). p-forms can be obtained from (antisymmetrized) wedge tensor product A of 1-forms. For 
instance, the wedge product between rj and the 1-form £ = z dx + 2 dy is the 2-form 



w 



rj A £ = 6 x 2 y 7 dx A dy — 5 y z dx A dz — 10 y dy A dz (3) 



Notice the absence of dx A dx = 0, since the wedge product is antisymmetric; we also use that dx A dz = —dz A dx. 
There are (™) linearly independent 2-forms dx 11 A dx v = dx^ ® dx u — dx u ® dx^ ((g) is the tensor product), which span 
the coordinate basis of 2-forms. Any 2-form can be written as 
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a = ^tt F dx 11 A dx v (4) 

where a^ v = —a ufi . In the former example it is zu xy = 6x 2 y 7 , zu yz = 0, vo xz — —byz, The factor 1/2! in Eq. (j3]) 
takes into account the fact that each independent element of the basis appears 2! times in the sum over /i and v. 

If a and j3 are 1-forms on a 3-dimensional manifold, then the components of the product a A /3 will look as the 
Cartesian components of the vector product in Euclidean space: 

a A j3 = (a x (3 y - a y /3 x ) dx A dy - {a x /3 z - a z /3 x ) dz A dx + {a y j3 z - a z /3 y ) dy A dz (5) 

Any p-form a and g-form /3 satisfy 

a A p = (-l) pq f3 A a (6) 

Thus 1-forms anti-commute but 2-forms commute, etc. 

The exterior derivative d is a nilpotent operator (cP = 0) . If d acts on a function / (0-form) the result is the 1-form 
df = {df /dx^) dx^. In general, if d acts on a p-form a then the result is a (p + l)-form da. Since d(dx fJl ) = 0, da is 
obtained by merely differentiating its components as exterior derivatives of functions: 



da = d ( «aa«/....1 a dxX A dx>1 A dx v .... (7) 

For instance 



p\ 



dr\ = -21 x 2 y 6 dx A dy + 5 dy A dz (8) 

If a is a p-form then 

d(a A 13) = da A 13+ {-If a A d/3 (9) 

A given p-form a is called closed if da — 0; it is called exact if it can be written as the exterior derivative of a 
(p — l)-form. Each exact form is closed; the inverse proposition is locally true, but its global validity depends on the 
topology of the space (see Poincare's lemmai^). 

Finally, we will introduce the Hodge star operator *. This operator changes p-forms into (n — p)-forms, and involves 
the components g^ v of the metric tensor present in the space-time interval. If a is a p-form whose components are 
a Ml ^ then 

* a„P+i = ^/|det(^)| £ Mp+1 Mp a pl ^ (10) 

where the indexes are raised with the inverse metric tensor g^ v : a Pl Mp = g^ lVl g^v v v a Vl Vp . e is the Levi- 

Civita symbol, which takes the value 1 (—1) for even (odd) permutations of the natural order of its indexes, and it is 
zero if there are indexes of equal value. The successive application of the Hodge star operator on a p-form a is££ 

* *a = (_l)P(»-P)+(™-<0/ 2 a (11) 

where a is the signature of the metric tensor (the difference between the numbers of positive and negative eigenvalues 
of the metric tensor). 



III. POTENTIAL AND FIELD AS DIFFERENTIAL FORMS ON A MANIFOLD 



The electromagnetic field is an exact 2-form F = dA, where the 1-form A is the potential. Given a set of non- 
necessarily Cartesian coordinates x, y, z together with the time t, one can express these forms in the coordinate basis 
{dt, dx, dy, dz}: 

A = A v dx" (12) 

If x, y, z are Cartesian coordinates, then the components A v will coincide with the scalar and vector potential: 
A v = (— <fi, A) (SI units) or A v — (— c<fi, A) (Gaussian units). Therefore 

F = dA = dA v A dx" = d„A u dx* A dx" (13) 

Since dx^ 1 A dx" = dx^ <S> dx" — dx" ® dx^, then the (antisymmetric) components of F are F^ — d^A v — d^A^. For 
Cartesian coordinates, it is F it = Ei in SI units (or cEi in Gaussian units), F yz — B x , F zx — B y , F xy = B z : 

F = E A dt + B x dy A dz + By dz A dx + B z dx A dy (14) 

where E ~ Ei dx 1 (SI units) is a 1-form. As Eq. (|14[) shows, in the 3-dimensional Euclidean space the wedge 
product between 1-forms works as the vector product, since dy A dz is the basis element for the a; -component of the 
pseudo- vector B, and so on (think the wedge product as a vector product; see Eq. ((U). 
Since F = dA is exact and d is nilpotent, then it results the identity 

dF = (15) 

which amounts those Maxwell equations used to define the potentials: V x E = —dH/dt, V ■ B = 0. 

The dynamical sourceless Maxwell equations, VxB = — c~ 2 dE/dt, V • E = 0, come from varying the electromag- 
netic action S[A] = — (4/i Q c) _1 J *F A F. Since F is a 2-form in a manifold of dimension n = 4 (the spacetime), then 
*F is a 2-form too. *F A F is a 4-form (a volume in spacetime). The resulting Eulcr- Lagrange equations are 

d*F = (16) 



IV. CYLINDRICAL WAVEGUIDES 



We are going to solve the sourceless Maxwell equations for the electromagnetic field in hollow cylindrical waveguides. 
So, let us call z the Cartesian coordinate along the waveguide, and x, y two properly chosen coordinates spanning the 
waveguide section: Cartesian coordinates for rectangular section, polar coordinates for circular section, etc. We will 
begin by working out stationary waves. Afterwards we will introduce the propagation along the waveguide by means 
of a Lorentz boost in the z -direction. Therefore we will start by proposing a solution independent of z. 



A. Stationary TM modes 

Let us consider a monochromatic potential having just a component along the waveguide: 

A = e int ip(x,y) dz (17) 

Function ip is defined in the waveguide section and has units of magnetic field times length. Thus the electromagnetic 
field is 

F = dA = e mt (ifl tp dt + dip) A dz 

= e m (ift ip dtAdz+ d x ip dx Adz + d y tp dy A dz) (18) 

The first term is an electric field along z and the other ones make up a magnetic field orthogonal to the z-axis, so lying 
on the waveguide section. Therefore the proposed solution is a TM mode. For the moment this solution does not 
propagate along the waveguide since the components do not depend on z (this feature will be introduced later). Field 
([T%1) is just a stationary wave bouncing between the boundaries. Function -0 is subjected to fulfill perfect conductor 
boundary conditions: the tangential electric field and the normal magnetic field must vanish on the boundary. Then 
ip must be zero on the boundary to cancel out the (pure) tangential electric field: 

1^\boundary (19) 
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(Dirichlet condition for the potential ip(x,y)). Since ip is constant on the boundary then the 1-form dtp = {dip/dx^) dx^ 
is normal to the boundary; so the magnetic part dip A dz is tangential to the boundary and perpendicular to the z 
-axis (just think the product as a vector product). Thus, the boundary condition for the electric field guarantees the 
boundary condition for the magnetic field too. 

The 2-form *F involves the metric tensor. If orthogonal coordinates are used in the waveguide section, then the 
space-time interval will read 

ds 2 = -c 2 dt 2 + g xx dx 2 + g yy dy 2 + dz 2 (20) 

where the metric di&g(g xx , g yy ) in the section depends on the choice of coordinates x, y. Then the determinant of the 
metric and the inverse metric tensor are 

I <iet(g^)\ = c 2 g xx g yy , g^ v = diag(-c~ 2 , l/g xx , l/g yy , 1) (21) 

To write *F we just need a few results: 

* (dt Adz) = -cT 1 ^Jg xx g yy dx A dy , 

*{dxf\dz) = -c^g xx g yy g xx dt A dy , (22) 

*(dy A dz) — c ^Jg xx g yy g vv dt A dx 

Thus, 

* F = c^g xx g yy e mt {-ic~ 2 Q ip dx A dy + ( g xx d x ip dy - g vv d y ip dx) A dt} (23) 



d*F = c 



y/g xx g yy e iQt (c~ 2 tt 2 ip + (2) A?/>) dt A dx A dy (24) 



Eq. (|16j) implies that function ip must be an eigenfunction ip mn of the two-dimensional Laplacian operator in the 
waveguide section, —c~ 2 fl 2 nn being its respective eigenvalue: 

1 £1 2 

(2) Aip mn = —= [d x (^g xx g yy g xx d x ip mn ) + d y (^Jg xx g yy g yv d y ip mn )] = ^ ip mn (25) 

\/ QxxQyy C 

Notice that Eq. (|2"3")l contains the two-dimensional 1-form ^ 2 \*dip) (the superscript "(2) " means that the Hodge 
star is applied in a n = 2 submanifold; see Appendix). As can be seen, the Laplacian in the waveguide section is 
( 2 K =( 2 )(— *d*d). Its eigenfunctions ip mn satisfying proper boundary conditions are identified by two discrete indexes 
m, n. Table U summarizes the results for monochromatic TM stationary modes. Table [TT] shows the solutions ip of 
Eqs. ([TO]). ([23)1 and their eigenvalues (allowed frequencies for stationary TM modes) for typical waveguide sections. 

TABLE I: TM stationary modes. 



F TM = e in mn t {iQ, mn ^ mn dt+ dip mn )Adz 
*F™ = e iUmnt (-ic _1 Sl m „ ^gxxg yy ipmn dx A dy + c (2 \*dijj mn ) A dt) 
ip and (2 \*dip) are denned in the waveguide section: 

^ ^Alp mn (x, y) = C £l mn Ipmn (x , y) Ipmn {boundary = 

metric in the section: d£ 2 — g xx (x,y) dx 2 + g yy (x,y) dy 2 



B. Stationary TE modes 



Remarkably F and *F are on an equal footing in Eqs. ([TBI) and (JTHJ) . This circumstance will allow to built the 
stationary TE modes by exchanging their roles. In Eq. (|14[) . this exchange amounts to the interchanging of E and 
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TABLE II: TM ip's and fi's for rectangular and circular sections. 

Rectangular section 

x, y are Cartesian coordinates, 0<x<a, 0<y<b, g xx = 1, g yy — 1 

■4) mn (x,y) = A m „sin (™x) sin (f-y) , m,n GN 

Circular section 

x,y are polar coordinates r,ip: < r < R, < (p < 2%, g rr = 1, g vv — r 2 

Ipmn (r, If) = Jm (imn T / R) (A mn COs(m<p) + B m n Sm(mVj)) 

fimn = Xm S c Xmn are zeros of Bessel functions J m 



B. So we define 

F TE = ^pTM (26) 

We will apply Eq. (fTTj) to obtain *F TE = * * F™ . In the space-time it is n — 4 and a — 2 (see Eq. ((20)) ); so 
** = Thus 

* F TB - -F™ (27) 

The field (|26|) can be ascribed to the potential 

pint „if2t„ 

^ = -TTTxra^ (ff ra ^ dx - 3 XX dy) = — (2 t*#) (28) 

In fact, by differentiating the middle expression in Eq. (|28[) one recognizes the appearance of the Laplacian defined 
in Eq. flS2J|: 

F T£ = dA TE = e int c ^/g xx g yy j ^ (2) A^ dxhdy+{ g xx d x ip dy - g vy d y i/j dx) A (ft 



c V9xxg yy { ^ (2) A^ cfe A dy +W(*dip) A <ft \ (29) 



This result coincides with field (|23[) when the pair f2,y> is chosen among the solutions £l m n,'4 , mn of the eigenvalue 
equation (|25j) . Although i 7,7 ^ solves both Eqs. (p~5|) and (fTBT) . the boundary condition should be consistently changed: 
the electric field should be normal to the boundary to accomplish the perfect conductor boundary condition. If 
n = n x dx + n y dy is a 1-form normal to the boundary, and ^ 2 \*dip) in Eq. (|29[) is a 1-form proportional to the electric 
field in the waveguide section, then one demands 

( 2 {*dlp) A Tl \boundar y = (30) 

(just think the product as a vector product). When written in vector language this requirement means 

n • Vtj) \boundary = (31) 

(Neumann boundary condition for the potential tp(x,y)). Table ILTI1 summarizes the monochromatic stationary TE 
modes. Table ITVl shows the functions tp for typical waveguide sections. 



V. PROPAGATING MODES 



The solutions studied in Section IV do not propagate energy along the waveguide. Since the Poynting vector is 
proportional to E x B, then there should exist E and B in the waveguide section to have energy propagating along 
the waveguide. However, field F™ in Eq. (|18l) has only a magnetic field in the waveguide section, and field F TE in 
Eq. (|2"§1) has only an electric field in the section. Thus, the Poynting vector in solutions (TT51) . (I2l)l) is orthogonal to the 
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TABLE III: TE stationary modes. 
F TE = e lf2m "* (-i c _1 n mn ^Jg xx g yy ipmn dx A dy + c ( 2 \*di} mn ) A dt 

*F TE = - e <f2 ">»< (ifimn V>mn rf* + #mn ) A C?Z 

^> and < - 2 {*d'0) are denned in the waveguide section: 

(2) A-i/w(s;,2/) = — C -2 fi 2 nn tpmn(x, y) n ■ jboundary = 

metric in the section: d£ 2 — g xx (x,y) dx 2 + g yy (x, y) dy 2 



TABLE IV: TE ?/>'s and Q's for rectangular and circular sections. 

Rectangular section 

x, y are Cartesian coordinates, 0<x<a, 0<y<b, g xx = 1 , g yy — 1 

ip mn (x,y) = A mn cos (™x) cos (^j/) 

Circular section 

a;,y are polar coordinates r,ip: < r < R, < ip < 2ir, g rr = 1, g vv — r 2 

1pmn(r, if) = J m {ymn T / R) (A mn COs(lTL(p) + B mn Sm(m(f)) 

£lmn = Vm S ° y-mn are zeros of the derivatives of Bessel functions J„ 



waveguide axis; so solutions (fT8|) . (|29|) are stationary waves where the energy bounces between the boundaries but 
does not propagate along the waveguide. Moreover, it is easy to prove that the time averaged Poynting vector vanishes 
in this case. This means that the fields (TT51) . (j2T))) are displayed in their proper frame. The stationary solutions (|18l) . 
(|29[) can be transformed into solutions that propagates energy along the waveguide by performing a Lorentz boost in 
the z-direction. To do this we use 

t = j(V)(t' - Vc~ 2 z'), dt = j(V)(dt' ~Vc- 2 dz'), dz = i{V){dz' -Vdt') (32) 

where j(V) = (1 - V 2 ^ 2 )- 1 / 2 . Thus 

dtAdz = 7(V) 2 (1 - V 2 C - 2 ) dt 1 A dz' = dt 1 A dz' (33) 

A. TM modes 

Eq. (|33p means that the longitudinal electric field remains invariant in Eq. (|18p . On the contrary the transverse 
term dip(x, y) A dz changes to 

dip(x, y) A dz = 1 {V) dtp(x, y) A [dz 1 - Vdt') (34) 

which amounts not only a change of the magnetic field in the waveguide section but the appearance of a transverse 
electric field. Of course this is nothing but the usual rules to transform electric and magnetic fields (see for instance 
Ref.[l6|). However the geometric language shows it in a quite elegant way. 
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B. TE modes 



In this case the Eq. (|33|) means that the longitudinal magnetic field remains invariant in Eq. (|29|) . Instead, the 
Lorentz boost changes the electric transverse sector of F TE to 

c {2 \*dtl) mn ) A eft = c 7(V) {2 {*d^ mn ) A {dt' ~ Vc~ 2 dz') 

= c 700 g yy d y ^ dx - g xx dy) A {dt' - Vc^dz 1 ) 

Therefore, not only the electric transverse field is changed by the boost, but a magnetic field appears in the waveguide 
section in the new frame. 

As a conclusion of this Section, in any frame differing from the proper frame where the solutions (fT5|) . (|29p were 
built, the propagating TM and TE modes display both electric and magnetic fields in the waveguide section and so 
they propagate energy along the waveguide. 



VI. TRANSMITTED POWER 



As already stated, the existence of both magnetic and electric transverse fields in the new frame produces a Poynting 
vector along the waveguide. Thus, in the new frame there is energy propagating in the waveguide. Velocity V is 
the velocity relative to the proper frame. In this sense V has the right to be called energy velocity, since no energy 
propagates along the waveguide in the original proper frame (there is just an energy flux orthogonal to the waveguide 
axis whose time averaging vanishes) . 

The time averaged energy flux along the waveguide, in the frame moving with velocity V relative to the proper 
frame, results from the t ' z' component of the electromagnetic energy-momentum tensor T^Jl 

Mo = F W F V p - - 9)iV F Xp F x P (35) 

Only real fields should be considered; so one has to average products of trigonometric functions. It is well known that 
< sin 2 (fM) >= i =< cos 2 (Sli) >, < sin(fit) cos(ftt) >= 0. Thus one obtains 

Ho < T Vz , >= g xx < F t , x F z , x > +g yy < F Vy F z , y >= ~^{V) 2 V |V^| 2 (36) 

(we remark that the Lorentz boost does not change the components of the metric tensor) . The result ([36]) is shared 
by TM and TE modes (although it is harder to get it for TE modes). 

c 2 T* 2 = —Tfi z i is the energy per unit of time and area going through the waveguide section. The transmitted 
power results from integrating this quantity in the section S. Since ip vanishes on the boundary of the section (its 
normal derivative vanishes at the boundary for TE modes), then the Green's first identity implies that^i 



• Vip dS = - J i> (2) Ai/; dS (37) 
(both for TM and for TE modes). Thus one uses Eq. (J2SJ) to obtain 

J | V</>| 2 dS = J c~ 2 n 2 *p 2 dS (38) 



Therefore the transmitted power is 



Pmn = I C 2 T tV dS = 7(V 9 )2y "- / An dS (39) 

J 2 yi c l J 

This result can be written in terms of the frequency u) and the wavenumber k z i . The transformation of the coordinate 
t implies that the wave phase becomes 

Qt = Sli(V)(t - Vrr 2 z') (40) 

Then 



lo = 7(V)n, k z , = 1 (V)VLVc- 2 (41) 
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which leads to the dispersion relation 

uj = yjc 2 k 2 z ,+Q 2 (42) 

(then Q mn is the cut-off frequency for each mode) . Thus the energy velocity V is written in terms of ui and k z i as 

c 2 k z i dui 
lo ok z i 

As could be expected, the energy velocity coincides with the group velocity duj/dk z i. 
Notice that wk z i equals the expression contained in the transmitted power. So 

Pmn = ^L[^ mndS (44) 
2/io J 



If the waveguide is filled with an homogeneous linear medium, then the Eq. f|44[) for the transmitted power can be 
used in the frame where the medium is at rest by replacing fi with the permeability /i (the constitutive relations are 
only valid in the media proper frames). 

To end the study of the energy transmission let us compute the time-averaged energy density c 2 T 4 * = c~ 2 T t 't' ■ 
The results to be obtained are common to TM and TE modes, but it will be easier to compute them for TM modes. 
We start from 

c 2 

fi < T t , t , >= g xx < F Vx F t , x > + g yy < F t , y F t , y > + < F Vz ,F t , z , > + — < F Xp F Xp > (45) 

The invariant < F\ p F Xp > can be computed with the components of the stationary field (it is invariant under a 
Lorentz boost!). The field of stationary TM modes has only three independent components. The real fields are 

Ft z — ^mn ^Pmn Sm(£7 mn £) 

F xz = d x tp mn cos(n mn t) (46) 

Fy Z dylpmn COs(H mn £) 

Therefore, the time-averaged scalar invariant < F\ p F Xp > is 

< F\ p F Xp >= 2 < F tz F tz + F XZ F XZ + F yz Fy z > (47) 
= -cr 2 VL 2 m J> 2 mn + <T (d x ^ m n) 2 + g yy {d v ^ mn ) 2 = -c- 2 n 2 m jj 2 mn + |VV„ < 2 



-'run. 



Then 

»o < t vv >= ^{vfv 2 ^ 2 + \v?i> 2 - \ n 2 ip 2 + ^|w| 2 = -n 2 ^ 2 + £ i±g |v^| 2 (48) 

c 2 

By performing the integral in the waveguide section we obtain the energy per unit of length: 

Umn = J c 2 T t ' t ' dS=^J j^L^n + \^ |V^ mn | 2 1 dS = J tf mn dS (49) 

where we have used the result (|38| . By comparing with Eq. (|44"]) one gets 

Pmn — V U mn (50) 

Once again we find the velocity V in the expected role of the energy velocity: in this case as the quotient of transmitted 
power and energy per unit of length as is usually defined. 6 Notice that T z 1 is the density of the z'-component of 
the momentum. Therefore, since the energy-momentum tensor is symmetric, the former expression can also be read 
as the momentum per unit of length equaled to the energy per unit of length over c 2 times the energy velocity 
(the usual relativistic relation for massive particles!). Even though the electromagnetic field is massless its energy 
in the waveguide behaves like the one of a massive field, as a consequence of the boundary conditions imposed by 
the waveguide. The same feature emerges in theories with compactified dimensions, which impose periodic boundary 
conditions to the (otherwise massless) fields.— 
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VII. APPENDIX 



Let (p,tp be two differentiable functions (0-forms). Then one has the following identity among volumes (n- forms), 

d(f> A *dip + (p d * dip — d((p * dip) (51) 

which can be integrated to become 



(d(p A *dip + (p d * dip) = (p * dip (52) 

S JdS 

(Stokes theorem^— has been used on the right side). The n-form d*dip is connected with = — *d*dip. According 
to Eq. (fTT|) it is * Atp = — d * dtp whenever ip is a 0-form and the space has signature a = n. Thus we get the Green's 
first identity 



(dep A *dip - <p * Atp) = <p * dip (53) 

S JdS 

To prove Eq. (|38l) we will apply Green's first identity to the case tj> = ip in the waveguide section, which is a 
2-dimensional manifold with coordinates x, y and signature a — 2. For any differentiable function (0-form) ip one has 

dip = d x ip dx + d y ip dy 

*dtp = y/g xx g yy ( g vy d y ijj dx - g xx d x ip dy) (54) 
dtp A *dtp = -^/g xx g yv (g xx (d x tp) 2 + g yy {d y tp) 2 ) dx A dy 

As stated, the 2-form dtp A *dip is a volume in the 2-dimensional waveguide section S. It contains the surface element 
dS = y/g xx g yy dx dy. Notice that the integral we are interested in Eq. (|38| is 

|VV>| 2 dS = - [ dtp A *dtP (55) 
Js 

So, we can compute it by using Eq. (|53[) for <p = ip. Since we are interested in functions ip that vanish on the boundary 
(TM modes), or their normal derivatives vanish on the boundary (TE modes) (i.e., *dip restricted to the boundary is 
null; see Eq. (|30)) ). then the right side of Eq. (|53"]) is zero in these cases. Besides, it is *1 = Jg X x9yy dx A dy, which 
will be used to compute the Hodge star of the 0-form Atp. Moreover, we are working with functions accomplishing 
the eigenvalue equation Atp — — fl 2 c~ 2 ip; then it is *Atp = — fl 2 c~ 2 ip ^/g xx g yy dx A dy. Thus the result is 

/ \Vip\ 2 dS =- dip A *dip = - / ip *AiP= / c~ 2 n 2 iP 2 y/g xx g yy dx A dy = / c~ 2 n 2 tp 2 dS (56) 
J Js Js Js J 



It should be emphasized that expressions like the ones in Eq. (|54[) only depend on the normalized basis of forms 
and vectors: 

e x = Vffer dx, e y = y/g^ dy (57) 
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